Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



O 

o 

(N 

■4— » 

o 
O 

m 



S3 
O 
O 

u 

Oh- 
o3 



13 
C 

o 

(N 
> 

m 
O 

o 

-I— > 
c3 



i 

C 

O 

o 



X 



Theoretical Study on Superconductivity in Boron-Doped Diamond 
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We consider superconductivity in boron (B) doped diamond using a simplified model for 
the valence band of diamond. We treat the effects of substitutional disorder of B ions by the 
coherent potential approximation (CPA) and those of the attractive force between holes by the 
ladder approximation under the assumption of instantaneous interaction with the Debye cutoff. 
We thereby calculate the quasiparticle life time, the evolution of the single-particle spectra due 
to doping, and the effect of disorder on the superconducting critical temperature T c . We in 
particular compare our results with those for supercell calculations to see the role of disorder, 
which turns out to be of crucial importance to T c . 
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1. Introduction 

Recent discovery of superconductivity in boron (B) 
doped diamond 1 ' 2 ) has renewed interest in physics of 
doped semiconductors. Inelastic x-ray scattering exper- 
iment 3 ^ indicates that the doped carriers are strongly 
coupled with bond-stretching phonons of u ~ 164 meV 
associated with the covalent bonds into which the holes 
are doped. The relatively high critical temperature of 
T c ~ 10K 4 - 5 ) is remarkable in view of the expected small 
density of states at the Fermi energy ep due to small car- 
rier density 6 - 1 as well as the existence of strong disorder 
intrinsic to the randomness in doped semiconductors. 7,8 ) 
In particular, the latter is deduced from the electrical 
resistivity measurement 9-11 ) and the line shape of quasi- 
particles observed by the angle-resolved photoemission 
spectroscopy (ARPES), 12 - 1 indicating that the life time of 
electrons r is very short with the characteristic parame- 
ter of epT < 1 (throughout this paper we set % = fee = 1) 
and violates the Ioffe-Regel criterion for coherent Bloch- 
like transport. 7, 8 ' n ) This small value of £ft implies that 
the Fermi surface is not well-defined in momentum space 
but is blurred by scattering (see Fig. 1); the supercon- 
ductivity in B-doped diamond may therefore be referred 
to as "superconductivity without Fermi surface". 13 ) 

The strong scattering potential of the impurity B ions 
may be responsible not only for such disorder-induced 
anomalous features but also for the basic electronic 
states. The energy level of a single B impurity in dia- 
mond is very deep, locating at 0.37 eV above the top of 
the valence band of diamond, 14,15 ) and, in the presence 
of such a strong scattering potential, the simple picture 
of the 'rigid-band shift' of the Fermi level due to doping 
(or the virtual crystal approximation) based on the first- 
principles LDA (local-density approximation) band cal- 
culations 16-18 ) with the adjusted chemical potential may 
not necessarily be justified. In fact, some spectroscopic 
experimental data 19-21 ) seem to suggest the presence of 
non-rigid-band-likc features reminiscent of the impurity 
states even in the metallic regime of B-doped diamond, 
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Fig. 1. Schematic representation of the Fermi surface in the two- 
dimensional analogue of the Brillouin zone for hole-doped dia- 
mond; (a) without disorder, a small Fermi surface is present, but 
(b) with disorder, the Fermi surface is smeared out by scattering. 



which has been focused upon by Baskaran. 7, 8 ) Hence, 
evolution of the impurity states due to doping should be 
accessed with more care. 

Another key feature is the observed very small car- 
rier concentration. A recent analysis of the NMR spec- 
tra 6 ) suggests that the real doping rate of carriers (or the 
number of carriers per site) is at most only ~ 1%, which 
differs very much from the nominal doping rates of ~ 5%, 
and that this is due to the formation of boron-hydrogen 
(BH) complex; i.e., H atoms inevitably introduced in the 
synthetic processes of the materials absorb the holes in- 
troduced by B impurities. This possibility has been sug- 
gested recently both theoretically and experimentally; an 
LDA calculation has shown that the B ions occupy the 
carbon sites substitutionally and the H ions sit on the 
interstitial positions between B ions and the neighbor- 
ing carbon ions, 22 ) and based on the assumption of the 
presence of these two types of B ions, the observed NMR 
spectra can be analyzed consistently. 6 ) We should there- 
fore reexamine the electronic states of B-doped diamond, 
in particular, in the regime of low carrier concentration 
of ~ 1%, starting from the dilute limit of B doping. 

In this paper, we will first study the electronic states 
of this highly disordered system, based on a simplified 
model for the valence band of diamond by use of the co- 
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herent potential approximation (CPA) for treating the ef- 
fects of the substitutional disorder of B ions. We thereby 
calculate the density of states, life time, and evolution 
of the single-particle spectra due to doping, the results 
of which are compared with available experimental data 
including ARPES spectra. 

We will also calculate the superconducting pairing sus- 
ceptibility based on the ladder approximation under the 
assumption of instantaneous interaction and estimate the 
doping dependence of the critical temperature T c . We 
thereby want to clarify the effects of disorder on the su- 
perconductivity in such unusual situations as described 
above. We will in particular compare our results with 
those for supercell calculations (i.e., the periodic arrange- 
ment of B ions) to see the role of disorder, which turns 
out to be of crucial importance. Here we do not pursue 
processes of the attractive interaction in detail, though 
we assume it is due to the electron-phonon interaction as- 
sociated with bond-stretching modes with very high fre- 
quency, as has been indicated both experimentally and 
theoretically. 3 ' 16 ~ 18 ' 23 ) Preliminary results of our work 
have been presented in ref. 24 ) 

This paper is organized as follows. In §2, we present 
our model and method of calculations. In particular, we 
derive the pairing susceptibility within the framework of 
CPA. In §3, we present our results of calculations for the 
density of states, life time, single-particle spectra, and su- 
perconducting critical temperature, which are compared 
with available experimental data in §4. Discussions are 
given in §5, and we summarize our work in §6. 

2. Method of Calculations 

2.1 Model 

In order to extract the effects of disorder on the elec- 
tronic states near the top of the valence bands of dia- 
mond, we assume, for the sake of simplicity, the tight- 
binding model based on the simple cubic structure. By 
this simplification, the degeneracy of three bands at the 
top of the valence bands of diamond is ignored. However, 
we believe that the essential features of the effects of dis- 
order can be taken into account in this simple model. 

Adding the random potential A at B sites and the on- 
site attractive interaction — V (V > 0), our model may 
be written by the Hamiltonian 

H = ^2(e k - c\ a c ka + A ^2 Uia ~ V Y n *T n a 

fcer i£B (7 i 

(1) 

with 

£k = t (cos k x a + cos k y a + cos k z a), (2) 

where rii a is the number operator at site i and spin a 
and c ka creates an electron with momentum k and spin 
(7, ji is the chemical potential. 

Hereafter, we assume t = 1 as the unit of energy un- 
less otherwise indicated; to reproduce the width of the 
valence band of diamond 22 eV, this choice corresponds 
to taking 3.7 eV as the unit of energy. The effective mass 
at the top of the valence band of diamond is roughly 
reproduced by this choice. 



0.30 - 




Fig. 2. Model density of states for the tight-binding band on the 
three-dimensional simple cubic lattice without disorder. Energy 
level of a single impurity is also shown. 

We determine the value of A so as to reproduce the 
energy level of B ions in the dilute limit, which is exper- 
imentally known to be located at 0.37 eV above the top 
of the valence band of pure diamond; 14 ' 15 > we obtain the 
value A = 2.30, which is used throughout the present 
work. The density of states of this model is shown in 
Fig. 2 with the energy level of a single B impurity. The 
sites with A are distributed randomly over the lattice 
sites of the simple cubic structure; this substitutional 
disorder is treated by CPA. We define the concentration 
of B ions as cb = N&/N where Nb is the number of B 
sites and N is the total number of sites in the system. 

2.2 Coherent Potential Approximation 

Let us assume the absence of the attractive interaction, 
i.e., V — 0, for the moment. In CPA, the Green's function 
is defined with the momentum-independent self-energy 
function T,(ie n ) as 

(fc, ie n ) = ie n - e k + [i - T,(ie n ) (3) 

where e k is the band structure without impurities eq. (2), 
H is chemical potential, and e n = (2n + l)n/(3, with 
(i = 1/T. The self-energy function in the presence of the 
substitutional disorder is given as the polynomial scries 
expansion 25-27 ) as shown in Fig. 3(a): 

oo 

£(ie„) = A Qs(cb) A 8 - 1 e^^ien) (4a) 

S=l 



1 - e(*e n )(A - E(ie n )) v ; 

where 

e(*£„) = ^^G(fc,ze„) (5) 
fc 

and Qs(cb) is the polynomial of c-q. In the following, we 
use abbreviations X = S(ie„) and O = Q(ie n ). 
The CPA equation for S 

9S 2 + (1 - A9)S - c B A = 0. (6) 

is solved numerically by use of eq. (2). 
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Fig. 3. (a) The self-energy S by CPA. Dotted and solid lines rep- 
resent the perturbation potential A and Green's function, re- 
spectively, and X represents Qs(cb). (b) Diagram for the pairing 
susceptibility El given in eq. (8). (c) Diagram for determining the 
vertex function T in eq. (10). (d) Diagram for the pairing sus- 
ceptibility x by the ladder approximation eq. (18). Wavy line 
represents the attractive interaction V. 



2.3 Pairing Susceptibility by CPA 

We consider the pairing (or Cooper-pair) susceptibility 
in the absence of V defined as 
r 

U(q,iu m ) = At e l ^ T (T T b\{T)b q ) (7) 
Jo 

with b q = J2k c -k-q/2,i c k- q /2,i and uj m = 2irm/p. 
Following the CPA theory for the magnetic susceptibil- 
ity, 28, 29 ) we find that eq. (7) is given by 



U(q,iu )m ) = Tj2J2 GG + T 



(8) 



n k 



with 



G = G(k,ie n ) (9a) 

G+ = G(—k + q, —ie n + iu m ) (9b) 

r = T(q,ie n ,iuj m ), (9c) 

which is illustrated in Fig. 3(b). The vertex function T 
is determined by the equation (see Fig. 3(c)): 

T = 1 + TA 



J2 & s Qs (e 5 ~ 2 + e s - 3 e+ + . . . + ef- 



(10) 



S=2 



with 



0+ = Q(-ie n + ioj m ) (Ha) 
A = A(q,ie n ,iu) m ) = ^-J2 GG +- ( llb ) 



Equation (10) is solved for F as 



r = i - a- 



with 



e-e +/ 

S + = E(-ie n + iu} m ). 



(12) 



(13) 



We thus obtain 

Il(q,iu; rn )=T^A(q ) (14) 



for the pairing susceptibility in the absence of attractive 
interactions. Since we have 

Q(-ie n ) - Q{is n ) 



A(0,fe„,0) = - 



2ie n + E(-ie n ) - E(ie„) 



(15a) 



F(0, ie n , 0) = 2?£ " + £( 7 £ " ) ' £(t£ " ) (15b) 

from eqs. (11) and (12), we see from eq. (14) that the 
pairing susceptibility for q = and iuj m — is given as 
follows; 



n(o,o ) = TV eH£ " ) ,- e(i£ " ) 



= ^/ d£/?(£) ^ tanh ( 



(3(e - n) 



(16a) 



(16b) 



where 



p(e) = -\ 5 Um Im 1 ^ G(fc, e - M + %$). (17) 
fc 

We thus find that the effects of disorder in the present 
approximation appear only through the density of states 
in disordered systems, p(s). 

We should note that the vertex correction T is essen- 
tial here; if we take into account the effects of disorder 
only from the self-energy via A, the value of the pairing 
susceptibility 11(0, 0) becomes very small (reflecting the 
effects of disorder too strongly), but the vertex correc- 
tion, which becomes very large around the Fermi energy, 
recovers the value of 11(0, 0), to result in the BCS form 
of the pairing susceptibility eq. (16). Thus, the effects of 
disorder, even though it is strong, can be renormalizcd 
solely into the density of states p(e). This cancellation 
is just as in the case of the calculation of the magnetic 
susceptibility. 28, 29) 

In the presence of the attractive interaction, V > 0, 
the pairing susceptibility at q = and iuj m — 0, x(T), is 
given in the ladder approximation by the process shown 
in Fig. 3(d), which results in 



X(T) 



11(0,0) 



(18) 



1-101(0,0) 

The critical temperature T c is determined from the con- 
dition V r II(0,0) = 1 at which x(T) diverges: i.e., 

'P{e-PY 



(19) 



Note that our result for 11(0, 0) can be considered as 
an extension of Anderson's theorem 30 ) which states that 
T c of an isotropic s-wave superconductor is not affected 
by nonmagnetic impurities if the density of states is 
unchanged by scattering. The extension is that, within 
CPA, even the strong disorder that deforms the density 
of states does not affect the validity of the BCS expres- 
sion eq. (16) for the pairing susceptibility. 

3. Results of Calculations 

3.1 Density of States and Self-energy 

Calculated results for the density of states are shown 
in Fig. 4. We find that the impurity states broaden and 
get weight continuously as the doping rate is increased. 
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Fig. 4. Density of states p(e) near the Fermi energy calculated 
by CPA. Local densities of states for carbon and boron sites are 
also shown. The vertical line indicates the Fermi energy ep in 
the absence of compensation. 



The broadening is due to the large imaginary part of the 
self-energy around the Fermi energy (see below). The im- 
purity states are separated from the top of the valence 
band up to the doping rate of cb — 0.2% (see Fig. 4(a)), 
but are merged into the top of the valence band with 
increasing doping. The weight of the impurity states cor- 
responds to the doping concentration. We assume the 
absence of compensation hereafter; i.e., the doping rate 
of holes is equal to Cb- Then, the Fermi energy ep at 
T = K (then fi = ep) is determined by 



f 

J — ( 



p(e)de = 1- ^. 



(20) 



Thus, the impurity states, if they are detached from the 
valence bands, are at half filling, as pointed out in ref.. 7 ) 



Fig. 5. (a) Diagram for defining the local density of states at B 
site, Pb(£)- Dotted and solid lines represent the perturbation 
potential A and Green's function, respectively, and X represents 
Qs(c B )- 



but the carbon component in the states is as large as the 
boron component. 



0.2 (a) c =0.002 




Fig. 6. Calculated real (upper panels) and imaginary (lower pan- 
els) parts of the self-energy E(e) by CPA. The vertical line indi- 
cates the Fermi energy in the absence of compensation. 



The real and imaginary parts of the self-energy S(e) 
are shown in Fig. 6, where we find that there is a large 
imaginary part — ImE(e) in the impurity states, which 
are responsible for the broadening of the density of states 
of the impurity states, as well as for the short life time 
of electrons (see below). 

3. 2 Life Time and Mean- Free- Path 




0.0 
0.00 



0.03 



0.06 
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The local density of states for B sites may be defined 
as in Fig. 5 and is given by 

PB (e) = --Im— V Q s (c B )A s e s (e — fi + iS) 

(21a) 

c B j m Q(e — fj, + iS) 

~ ~~ m 1 - (A - S(e - Ji + i5))0(e - Ji + id) 

(21b) 

with 5 — > + . The local density of states for carbon sites 
may be defined as pc{z) = p{s) — Pb{s)- The calculated 
results are shown also in Fig. 4, where we find that the 
boron component is concentrated on the impurity states 



Fig. 7. Calculated life time E-pr (•) and mean-free-path kp£ (o) 
as a function of cb- 

The life time r at the Fermi energy e F may be defined 

as 



2ImE(e F + i(5) v ' 

with 8 — > + , the value of which is determined directly 
from the imaginary part of the self-energy calculated at 
£p in the presence of disorder. By use of this r, the mean- 
free-path I — v-pT, EpT, and k-pi may be defined with 
«f, fcp, and Ep being the Fermi velocity, Fermi wavevec- 
tor, and Fermi energy, respectively, corresponding to the 
same carrier number but in the absence of disorder. 
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The calculated results for E-pT and kp£ are shown in 
Fig. 7 as a function of Cb, where we find that r and £ are 
indeed very short, irrespective of the doping rates, which 
comes from the large imaginary part of the self-energy 
at the Fermi energy (see Fig. 6). At cb = 0.01, e.g., we 
have the values EpT ~ 0.48 and kpi ~ 0.55. 
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with /(e) = l/(eP( e ~W + 1). Here, we assume the De- 
bye cutoff of lud = 2300 K corresponding to the bond- 
stretching mode; i.e., the energy integral of eq. (16) is 
from e = fx — wd to /i + LUp>. The calculated results for T c 
are shown in Fig. 9 as a function of cb for various values 
of the attractive interaction V. We find that the value of 
T c increases rapidly with doping rate cb and the attrac- 
tive interaction V. If the energy of the bond-stretching 
mode of diamond lo — 164 meV is used as the cutoff en- 
ergy in the integral of eq. (16), the value of T c is found 
to decrease by ^30% but the general trend seen in Fig. 9 
does not change. 
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Fig. 9. Superconducting critical temperature T c as a function of 
cb calculated by CPA. 



Fig. 8. (Color online) Calculated results for the single-particle 
spectra p(fc,e) at cb = 0.01 (upper panel) and cb = 0.03 (lower 
panel). The artificial broadening of S = 0.13 is included in the 
spectra for the illustrative purpose. 



3.3 Single- Particle Spectra 

The calculated single-particle spectra are shown in 
Fig. 8, where we illustrate the quantity 

p(k,e) = --ImG(k,e + iS), (23) 

by which we can envisage the broadening of the band 
structure due to disorder in the Brillouin zone as well 
as the location of spectral weight due to the impurity 
states. We find that the impurity states are clearly seen 
just above the top of the valence band. We note that the 
impurity states are fairly localized in momentum space 
around the T point even at very low doping rates. This re- 
sult will be compared later with the results of the ARPES 
experiment. 12 ^ 

3.4 Superconducting Critical Temperature 

The superconducting critical temperature T c is calcu- 
lated from eq. (19) together with the chemical potential 
y, determined from 

£b 
2 



de/(e)p(e) 



1 



(24) 



Our results may be compared with results of the first- 
principles LDA calculations based on either virtual crys- 
tal approximation 16, 17 ' or supercell method. 17,18 -* We 
find that our value of the density of states at the Fermi 
energy p(sp) is about a factor 5 smaller than the values 
of the LDA calculations 16 - 1 if we compare at doping rate 
of cb = 0.03. The causes of this difference are the fol- 
lowing; the absence of the three-fold degeneracy on one 
hand and the effects of disorder on the other hand in 
the present calculations. Actually, for the latter, we have 
the value p(e F ) = 0.027 (0.046) at, e.g., c B = 0.01 (0.03) 
in our CPA calculation, while without disorder in the 
same model, we have p(sf) = 0.035 (0.053) at cb = 0.01 
(0.03). 

The strength of the electron-phonon coupling constant 
A has also been estimated by the LDA calculations, 16 ~ 18 ) 
where it was found that the contribution stems predom- 
inantly from the bond-stretching phonon mode and that 
the values are A ~ 0.3 - 0.55 at c B ^ 0.02 - 0.03. The 
value of T c has thereby been estimated based on the BCS 
or McMillan's formula 32 ' and found to be in quantita- 
tive agreement with experiment. 16 ~ 18 ' Then, from the 
comparison with the results of the LDA calculations, the 
strength of the attractive interaction V in our theory 
may be estimated; i.e., by using the relation 32 ' 

p(sf)V = (25) 
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with our calculated value of p(sf) = 0.027 at cb = 0.01 
and a standard value of the Coulomb pseudopotential 
parameter /i* ~ 0.1, and by assuming a typical value 
of the electron-phonon coupling constant A ~ 0.31 at 
cb = 0.01, we can estimate the value of the attractive 
interaction to be V ~ 6. 

4. Comparisons with Experiment 

4-1 Density of States 

The calculated results for the density of states (see 
Fig. 4) may be compared with results of a number of 
spectroscopy experiments such as the photoemission and 
absorption spectra 12, 19, 20,31 ) and optical conductivity 
spectra. 21 ) Although the effects of the matrix elements 
relevant to each experiment is not taken into account in 
the present calculations, the comparison can offer impor- 
tant information on the impurity states of the hole-doped 
diamond. 

We first point out that our calculated results for the 
density of states (see Fig. 4) indicate that the impurity 
states evolve continuously by doping from the low-doping 
limit, where it has been established 14,15 -* that the triply- 
degenerate impurity level exists for a single B impurity in 
otherwise pure diamond. This is in contrast to the claim 
of the rigid-band shift of the Fermi energy suggested in 
the first-principles LDA band calculations by the super- 
cell method, 17,18 '' which cannot reproduce the separated 
impurity level at low doping limit of diamond. 22 ' 

The presence or absence of the impurity states in the 
metallic regime of B-doped diamond should in principle 
be clarified by experiment. First, the angle- integrated 
photoemission spectroscopy experiment 3 1 ' has not so 
far provided evidences for the presence of the impurity 
states. On the other hand, the X-ray absorption (XAS) 
and emission (XES) spectroscopy experiments 19, 20 ) seem 
to suggest the presence of the impurity states 8 - 1 although 
not yet fully conclusive. The optical conductivity experi- 
ment 21 -* also seems to support the presence of the impu- 
rity states even in the metallic regime, where it has been 
found that, by doping, there appear the broad spectral 
features at ~ 0.37 eV that may correspond to the impu- 
rity states while the Drude weight is hardly visible. 

Further experimental studies are desired to clarify the 
nature of the impurity states in B-doped diamond. 

4-2 Life Time and Mean- Free- Path 

The calculated results for the life time and mean-free- 
path are shown in Fig. 7. At cb = 0.01, e.g., we find 
the values E-pr = 0.48 and kpt = 0.55. These values of 
E-pT and kp£ are close to the experimental values Epr — 
0.27 and k F £ = 0.32, 0.48, and 1.0 for the three bands 
obtained from the line shape of quasiparticles by ARPES 
experiment for the highest T c sample. 12 ) The analyses 
of the measured resistivity 2, 9 ~ n ) lead to similar small 
values of E$t and kp£ as well. 

4-3 Angle-Resolved Photoemission Spectra 

The calculated results for the single-particle spectra in 
Fig. 8 may be compared with the experimental ARPES 
spectra given in Fig. 2 of ref.. 12 ) To do this, we recal- 
culate the spectra by convoluting calculated results with 




Fig. 10. (Color online) Calculated results for the single-particle 
spectra p(fc,e) at doping rates of (a) cb = 0.2%, (b) cp, = 0.8%, 
(c) Cb = 1.0%, and (d) cb = 1.5%. The finite resolutions in 
both energy (250meV) and momentum (±0.026au _1 ) in exper- 
iment 12 ' are taken into account in the calculation. The spectra 
below the Fermi energy (shown by horizontal line) may be com- 
pared with the experimental ARPES data obtained by Yokoya 
et al. 12 ' The unit of energy is given in eV for comparison with 
experiment. Horizontal axes are in units of 1/a. 



the experimental resolutions in both energy (250 meV) 
and momentum (±0,026 au _1 ) in view of the actual ex- 
periment. 12 -' The results near the Fermi energy are shown 
in Fig. 10, where the spectra both below (which can be 
compared with experiment) and above the Fermi energy 
are presented. We should stress again that the three-fold 
degeneracy of the top of the valence band is ignored in 
our calculations. We find however that the agreement be- 
tween theory and experiment (below the Fermi energy) 
is obviously reasonable; i.e., even in the presence of the 
impurity states, the spectra of the valence band indeed 
look just as if the rigid-band shift of the Fermi energy 
were operating. But this apparent agreement does not 
imply that the impurity states are not present. This is in 
particular the case at low-doping regions; at cb = 0.2%, 
the impurity states are not clearly identified in the calcu- 
lated angle-resolved spectra (see Fig. 10(a)), but in the 
calculated density of states (see Fig. 4(a)), we find that 
the impurity states certainly exist, which are separated 
from the top of the valence band. Thus, up to the dop- 
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ing rates of cb — 1%, the shape of the angle-resolved 
spectra below the Fermi energy is not largely modified 
by the presence of the impurity states. We therefore con- 
sider that the claim 12 ) of the absence of the impurity 
states made from their ARPES data may not be justi- 
fied, at least with resolutions in both momentum and 
energy used in the actual experiments. Further exper- 
iments with much higher resolutions would be able to 
resolve the issue. 




Fig. 11. Comparison between the calculated critical temperature 
T c as a function of cb and the experimental data (open circles 
and broad curve) estimated by Mukuda et al. e ' from the analysis 
of the NMR spectra. T c is given in units of K for comparison with 
experiment. The values of V used for calculations are given in 
units of eV here. 



4-4 Superconducting Critical Temperature 

The calculated results for the doping dependence of 
T c are compared with experiment in Fig. 11. We find 
that if we assume an appropriate value of V, e.g., V ~ 6 
(which corresponds to ~20 eV), our result of T c ~ 7 K at 
cb = 1% is in reasonable agreement with results of the 
LDA calculations 16-18 ) and therefore with experiment. 
In particular, we find that the rapid increase in T c with 
increasing doping rate, which is observed in experiment, 
is reproduced fairly well. 

If seen more quantitatively, however, one may notice 
that the calculated concave dependence of T c on cb (see 
Fig. 11) is somehow not consistent with experiments 
where the nearly linear dependence on cb in the loga- 
rithmic scale is seen. We consider that the discrepancy 
may partly be due to the insufficiency (or mean-field na- 
ture) of CPA. In CPA, both the width and height of the 
density of states of the impurity states vary as ~ ^/cb 
(hence, the area varies as ~ cb), whereas in reality, there 
is a tendency 33-35 ) that the width varies as ~ cb and 
the height is almost kept constant ~ 1. If we take into 
account this tendency in our calculations, the agreement 
between theory and experiment will be improved. 

5. Discussions 

In this section, we want to discuss some issues and 
possibilities. 

First, let us consider possible improvements of our cal- 
culations presented in this paper. We should first of all 
point out that our calculated values of T c are not of 
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sufficient quantitative significance since the calculations 
are based on the single-band model of the simple-cubic 
lattice and the instantaneous interaction for the attrac- 
tive force. For more quantitative estimation of T c , we 
may need the following improvements: (i) The nearly- 
degenerate three-band structure around the top of the 
valence band of diamond in the actual lattice structure 
should be included, (ii) The retardation effects for the 
clcctron-phonon coupling should be taken into account, 
(iii) The dynamical aspect of the Coulomb interaction 
between holes should be important; due to the effects of 
diffusion by randomness, the Coulomb pseudopotential 
parameter fj,* should be effectively larger. 36 ' 37 ) (iv) Ap- 
proximations beyond CPA will be desirable as pointed 
out in §4.4. These subjects deserve further studies. 



(a)c B = 0.008 (b)c B = 0.016 




(c)c =0.008 (d)c =0.016 




E £ 



Fig. 12. Densities of states near the Fermi energy obtained by 
CPA (upper panels) and those calculated under the assumption 
of the periodic arrangement of B ions (lower panels) . The 5x5x5 
and 4x4x4 supercells are assumed in (c) and (d), respec- 
tively. Partial densities of states are also shown where '6 n.n.C 
means the partial density of states summed over the six nearest- 
neighbor carbon sites around a boron ion. Vertical lines indicate 
the Fermi energy ep. 

Next, let us consider what happens to T c if there were 
no effects of disorder of B ions. To see this, we carry 
out the supcrccll calculations for the same set of param- 
eter values used in the CPA calculations. The B ions 
are arranged periodically by using 6x6x6, 5x5x5, 
4x4x4, and 3x3x3 supercells of the simple-cubic lat- 
tice, which correspond to 0.46%, 0.80%, 1.56%, and 3.7% 
doping rates, respectively. The calculated results for the 
density of states are shown in Fig. 12, which are com- 
pared with the results of CPA at the same doping rates. 
We find that a rather narrow density of states for the im- 
purity states appears above the top of the valence band, 
which form bands (impurity band), and that the density 
of states at the Fermi energy is much higher than the re- 
sult of the CPA calculation. The separation between the 
impurity band and the top of the valence band is ~ 1 
eV, which is significantly larger than the energy of the 
impurity level in the dilute limit, 0.37 eV, due to level 
repulsions. 
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Fig. 13. Superconducting critical temperature T c obtained under 
the assumption of the periodic arrangement of B ions. 



We then calculate the values of T c for these super- 
cell systems by using eq. (19). The results are shown in 
Fig. 13. We find that the obtained values of T c are nearly 
an order of magnitude higher than the results of the CPA 
calculations for the same attractive strength (see Fig. 9). 
Thus, we may conclude that, if impurity B ions could be 
arranged periodically, one would be able to have the crit- 
ical temperature of an order of ~100 K. This value of T c 
exceeds T c ~ 30 — 40 K of what is called the theoretical 
limit of the BCS theory for conventional superconduc- 
tors; this is not surprising because the Debye frequency 
of diamond wd = 2300 K is exceptionally large, together 
with the presence of the strong electron-phonon coupling 
because of the hole doping into the covalent bonds. 13 ^ In 
other words, we can say that, in actual B-doped dia- 
mond, the effect of disorder indeed strongly suppresses 
the value of T c . Note the 'convex' dependence of T c on 
cb in Fig. 13, which is in contrast to that of the CPA cal- 
culation where the dependence is 'concave' (see Fig. 9); 
i.e., the T c in the supercell calculations decreases very 
rapidly with decreasing cb- This comes from the narrow 
density of states of the impurity band, where the width 
of the impurity band becomes less than the Debye cutoff 
and the integral in eq. (19) becomes very small. 

Together with reduction of the effects of disorder to 
raise T c , the increase of hole concentration is obviously 
another way. In the actual experiment, one may try to 
find out way to remove or reduce the H atoms introduced 
in the synthetic processes of the materials; i.e., the re- 
duction of passivation rate. 

Finally, we want to point out the possibility of the oc- 
currence of superconductivity in similar covalent-bonded 
materials. Based on the present study, we may suggest 
that the doping of carriers in boron-nitride (BN) and 
silicon-carbide (SiC) may be quite interesting; these ma- 
terials with the cubic (zinc-blende) structure not only 
have the valence-band structure with a large band gap 
similar to diamond but also have the bond-stretching 
phonon mode with very high frequency. 38 ~ 44 ) Then, we 
may expect that the doping of enough carriers would lead 
to superconductivity with appreciable T c . 



6. Conclusions 

We have studied the superconductivity in B-doped di- 
amond by use of CPA for treating the substitutional 
disorder of B ions and by treating the attractive force 
between holes by the ladder approximation under the 
assumption of instantaneous interaction. The main mes- 
sages obtained from our studies are the following: 

(i) The effects of disorder arc essential for understand- 
ing the electronic state of B-doped diamond. The broad- 
ening of the density of states and single-particle spectra 
near the Fermi energy, as well as the very short life-time 
and mean-free-path, are due to the effects of disorder, 
i.e., the large imaginary part of the self-energy induced 
by the disorder of B ions. 

(ii) The pairing susceptibility derived by CPA takes 
the form of the BCS theory, where the effects of disorder 
are renormalized solely into the density of states, even in 
the presence of strong disorder that deforms the shape of 
the density of states. The vertex correction is found to be 
essential here: i.e., the susceptibility becomes too small 
if we take into account only the effect of the self-energy, 
but it is recovered by the vertex correction, leading to 
the BCS expression for the pairing susceptibility. 

(iii) The superconducting critical temperature T c can 
be raised substantially by reducing the effects of disorder 
in doping processes, e.g., by controled spatial arrange- 
ment of B ions. 

We believe that the doped semiconductors should have 
capabilities for producing a new class of superconducting 
materials and therefore we hope that our work presented 
here will encourage further experimental studies of su- 
perconductivity induced by doping into covalent-bonded 
semiconductors. 
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